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COMPUTER PREDICTIONS OF PHOTOCHEMICAL OXIDANT LEVELS 


FOR INITIAL PRECURSOR CONCENTRATIONS CHARACTERISTIC 

OF SOUTHEASTERN VIRGINIA 

N. T. Wakelyn and Allen G. McLain 
Langley Research Center 


SUMMARY 

A computer study has been performed with a photochemical box model, 
using a contemporary chemical mechanism and procedure, and a range of initial 
input pollutant concentrations thought to encompass those characteristic of 
the Southeastern Virginia region before a photochemical oxidant episode. The 
model predictions are consistent with the expectation of high summer after- 
noon ozone levels when initial nonmethane hydrocarbon (NMHC) levels are m 
the range 0 . 30 - 0. 40 ppmC and N0 X levels are in the range 0.02 - 0.05 ppm. 
Calculations made with a Lagrangian model, for one of the previously calcu- 
lated cases, which had produced intermediate afternoon ozone levels, suggest 
that urban source additions of NMHC and N0 X exacerbate the photochemical 
oxidant condition. 


INTRODUCTION 

The Southeastern Virginia region regularly experiences high photochemical 
oxidant concentrations (refs. 1 and 2) during the summer months. This 
condition is characterized by a general southwesterly flow of air into the 
region and bright "blue sky" days occurring m conjunction with a high- 
pressure cell. The constancy of summer weather presents an opportunity for 
the test of airborne and satellite instrumentation m a large and predictable 
area containing rural, urban, and water regions. The Southeastern Virginia 
Urban Plume Study of the summer of 1977 (ref. 3) was instituted to utilize the 
area as a test bed for remote sensors currently under development by the NASA 
and its contractors. A variety of ozone levels are measurable m the region 
from early morning levels as low as about 0.030 ppm to late afternoon levels 
which can reach, from time-to-time , almost 0.200 ppm (refs. 2 and 3). 

In order to fully utilize the potential of the data from such a study, 
photochemical models of the region must be employed to test the reasonableness 
and consistency of the various inputs to the data base. As a beginning effort 
in the development and test of a photochemical model, a finite rate chemical 
kinetics package was modified to vary the rate coefficients of the unimolec- 
ular photolytic decompositions in a diurnal manner. The package is capable of 
handling unimolecular , bimolecular, and termolecular chemistry and is provided 
with a stiff ordinary differential equation solver (ref. 4). Provision for 
time-varying source inputs of pollutant species, mixed-layer rise dilution, 
and restart and off-line plotting capability was also written into the 
computer code. 







A computer exercise was performed using a contemporary, chamber-verified 
chemical mechanism (refs. 5 and 6), a box- type model similar to that 
described in reference 7, and a test procedure (ref. 8) of merit. Input for 
the exercise spanned that range of common pollutant species concentrations 
generally thought to be characteristic of the Southeastern Virginia area 
during the morning preceding an afternoon photochemical oxidant incident. 
Preliminary calculations were also made with a Lagrangian (moving air parcel, 
ref. 8) model to determine the effect of a range of consistent source inputs 
on oxidant levels for the region. The purpose was to test the compatibility 
of the mechanism, model, and procedure and to obtain predictions and trends 
which could be compared to data obtained from future field experiments. A 
family of curves was generated of predicted ozone levels as a function of 
initial N0 X and nonmethane hydrocarbon (NMHC) input as well as source (time 
varying) inputs of these substances, with diurnal variation of the photolytic 
rates. The present paper is concerned with the results of these calculations. 


THE CHEMICAL MECHANISM 

The chemical kinetic mechanism used in this computer simulation consists 
of photochemical decompositions, inorganic rearrangements, organic oxidation 
reactions (with "lumping" of organics into paraffins and olefins), and 
reactions describing the rearrangements of oxidized organics. A rather 
complete sensitivity analysis is presented in reference 6. The following 
unimolecular photochemical decompositions, with the rate coefficients 
calculated by the method of reference 9 and varied every 5 minutes over the 
course of the simulation, were used: 


k,* min 


1 . 

N0 2 NO + 0 

2.036E-1 

2. 

hno 2 


OH + NO 

1.182E-2 

3. 

hno 3 


N0 g + OH 

2.233E-6 

L. 

H 2°2 


OH + OH 

U.520E-U 

5. 

RCHO 

->• 

stable products 

1 . 59^E-3 

6. 

RCHO 


1/2 R0 2 + 3/2 H0 2 

1.59UE-3 


*0800 l.d.t. (local daylight time), Hampton, Virginia, latitude 36.85 
Longitude 76.35, August k, 1977, "blue sky." 


The photolytic rate coefficients and their diurnal variation over a 10-hour 
period are shown in figure 1. 
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The inorganic rearrangements (from ref. 6) used are as follows: 


k, ppm mm (x = 0,1 ) 


7. 

o + (o 2 ) -»■ o 3 

3.30E+6 

8. 

0 3 + NO -► N0 2 + 0 2 

2.10E+1 

9. 

0 3 + no 2 -> no 3 + o 2 

4.60E-2 

10. 

N0 3 + NO + 2 N0 2 

1.50E+4 

11. 

N0 2 + N0 3 + (HgO) -► 2 HN0 3 

4. 50E+0 

12. 

0 + NO -* N0 2 

2.50E+3 

13. 

0 + N0 2 ■> NO + 0 2 

1.40E+4 

14. 

o + no 2 •> no 3 

3.00E+3 

15. 

NO + N0 2 + (H O) -»■ 2 HN0 2 

1.50E-2 

16. 

2 HN0 2 •> NO + N0 2 + H 2 0 

1.10E+0 

17. 

OH + N0 2 -> HN0 3 

6.00E+3 

18. 

OH + NO + HN0 2 

4.80E+3 

19- 

H0 2 + NO -»■ OH + N0 2 

7.00E+2 


The oxidation of hydrocarbons to produce aldehydes and alkyl and acyl 
peroxide radicals is described by the following reaction set for olefins 
(ref. 6) 


k , ppm min 


20. OLEF + OH + RCHO + R0 2 2.50E+4 

21. OLEF + 0 -> R0 2 + 1/2 hO^ + 1/2 RC0 3 6.80E+3 

22. OLEF + 0 3 -+■ RCHO + RC0 3 + OH 1.60E-2 

and, for paraffins (ref. 5)‘ 

23. PARA + OH -*• R0 2 + H O 3.80E+3 

24. PARA + 0 ■+■ R0 2 + OH 6.50E+1 
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The hydrocarbon mechanism was described by the referenced workers as smog 
chamber verified and specialized for automobile exhaust. As one would 
expect, the rate coefficients for reactions involving olefins are greater 
than those involving paraffins (olefins > paraffins) and OH > 0 > 0^ . 

The rearrangement of oxidized organics to produce, in part, PAN and 
alkyl nitrite and nitrate is described by the following set (ref. 6) 

— X — ~| 

k, ppm mm (x = 0,l) 


25. 

RCHO + OH + 1/2 RC0 3 + 1/2 HC> 2 

2.30E+4 

26. 

R0 2 + NO -*• R0 + N0 2 

3.00E+3 

27. 

RC0 3 + NO -* R0 2 + N0 2 

1.50E+3 

28. 

RCC> 3 + N0 2 -* PAN 

5.00E+2 

29. 

R0 + (02 ) -* RCHO + H0 2 

5.00E+3 

30. 

R0 + N0 2 -> R0N0 2 

4 . 90E+2 

31. 

R0 + NO -+ R0N0 

2. 50E+2 

32. 

2 H0 2 -*■ H 2 0 2 + 0 2 

5.30E+3 

33. 

R0 2 + H0 2 R0 + OH + 0 g 

1.00E+2 

34. 

2 R0 2 -> 2 R0 + 0 2 

1.00E+2 


MODELS AND TEST PROCEDURE 

There are tropospheric models in current use or under development of 
varying complexity in their description of atmospheric diffusion and 
chemistry. For example, see reference 10. In the event of inadequate 
meteorological data (a common occurrence) or for the purpose of emphasizing 
the chemistry part of the problem with a minimum of mathematical difficulty, 
a box (ref. 7), multi-box (ref. ll), or traveling air parcel (ref. 8) model 
can be used as the modeling tool. Since the purpose of the current computer 
study is to test the compatibility of a contemporary chemical scheme and 
test procedure with a range of initial input data for a region of interest, 
and to find trends which could help plan and interpret future data gathering 
studies, the box concept was used for the bulk of the work with some 
calculations performed with a traveling air parcel model. 

The procedure used in the computer simulation, essentially that of 
reference 8, is as follows 

a. Nonmethane hydrocarbons were taken to be a mix of olefin (propylene) 
and paraffin (n-butane) in such proportion as to simulate the reactivity of 
automobile emissions 25 percent propylene (ppmC) + 75 percent n-butane 
( ppmC ) . 
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b. Aldehydes were assumed to be 5 percent of the initial NMHC levels 
( m ppmC ) . 

c. All computer simulations were 10 hours long, beginning at 0800 l.d.t. 

d. Photolytic rate coefficients were varied every 5 minutes in accordance 
with diurnal variation in sunlight from 0800 - 1800 l.d.t. Mixed-layer dilutions 
method of reference 9 (see fig. l). 

e. An atmospheric dilution rate, due to mixed layer rise, of 3 percent 
per hour was assumed for the period 0800 - 1500 l.d.t. Mixed-layer dilutions 
for S.E. Virginia are generally substantially greater than this. 

f. Following the suggestion of reference 12. 

{ NO = 0.25 NO 

^ x 

NO = 0. 75 NO 

x 

Calculations were also performed without atmospheric dilution and with a 
different ratio of N02 and NO for N0 X in order to illustrate the magnitude 
of these effects. (See appendices A and B. ) 


Photochemical Box Model 

The phot oc hemic ally reacting volume under consideration in the box model 
is conceived as being rectangular in form, length not necessarily equal to 
width, and with a varying height to account for mixed layer rise. Since no 
provision for ground scavenging of ozone is provided m the present form of 
the model, the bottom of the box is at a position some 10 meters above ground 
level. The box is oriented such that the front face is normal to the wind 
direction, the initial concentration of pollutants and precursors at the 
front face is characteristic of the box as a whole, and any calculated 
concentration changes are homogeneously applied over the entire box. 
Photochemistry is calculated by marching forward in time from 0800 l.d.t. while 
changing the photolytic rate coefficients and outputting results every 5 
minutes . 

The box model used m this study does not consider sources within the 
modeled region and thus, by implication, only applies to over water and rural 
areas north of Norfolk. The traveling air-parcel model to be described later 
does not suffer from this inadequacy. Aside from the fact that current 
knowledge of local sources of N0 X and NMHC is less than satisfactory, there 
is need to know what ozone levels can result from the initial input alone. 

This baseline value of what levels of oxidant are to be expected after several 
hours of photochemical reaction, irrespective of local precursor contribution, 
is the concern of the major portion of this study. 
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Traveling Air Parcel Model 

To determine the effect of source loadings on eventual ozone levels, 
several calculations were performed with what, in essence, is a Lagrangian 
(traveling air parcel) model. The same procedure and photochemical 
mechanism was used for both the present model and the previously discussed 
photochemical box model. A parcel of air, with initial N0 X and NMHC 
concentrations of one of the intermediate (in ozone generation) cases 
previously calculated with the photochemical box model, enters the Norfolk 
area from the southwest, is in the Norfolk area and subject to source 
additions of N0 X and NMHC for a 2-hour period (0800 - 1000 l.'d.t.), and then 
passes out over the water where sources are "turned off" and where it 
pnoto chemically ages for the rest of the simulation. The 2-hour period 
is consistent with the urban area size and prevailing wind speed. 
Proportionate source contributions, with instant homogeneity assumed, were 
made every 5 minutes of the Lagrangian simulation. In keeping with the 
suggestion of reference 7, the following ratio is used for primary 
emissions of oxides of nitrogen: NO /NO = 0.10 and N0/N0 = 0.90. 

c. X X 

RESULTS AND DISCUSSION 
Photochemical Box Model 

The range of pollutant species concentrations thought to be 
cnaractenstic of the Southeastern Virginia region on a summer's morning 
before an afternoon photochemical oxidant incident is depicted in table I. 

A NMHC N0 X ratio of from 5 to 20 is considered characteristic of an urban 
area (ref. 12). 


X ' S \NMHC 

NO 

x \ 

0.10 

0.20 

0.30 

0 Uo 

0.050 

2.0 

u.o 

6.0 

8 0 

0.01+0 

2.5 

5.0 

7.5 

10.0 

0.030 

3.3 

6.7 

10.0 

13.3 

0.020 

5.0 

10.0 

15.0 

20 0 


Table I.- Range of initial NMHC (ppmC) and N0 X (ppm) 
levels used m simulation and their 
associated IJITIC . II0 x ratios 


Computer simulations were made for the l6 combinations of N0 X and 
nonmethane Hydrocarbons , as initial conditions, with II0 x varying from 
0.020 - 0.050 ppm and NMHC varying from 0 10 to 0.1+0 ppmC . The 
cnaractenstic incoming ozone is typically from 0.030 - 0.060 ppm (ref. 2). 
The simulations were based upon an initial ozone value of 0.050 ppm. 
Variations of ozone concentration with time for four initial NO incuts are 
snown m figure 2 (a, b, c, d), with each figure representing a aifferent 
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initial NMHC concentration. What appears to he a dip in the ozone 
concentration at 1500 hours for all the curves is the point in time at which 
mixed layer dilutions are no longer taken. For all of these simulations, 
there is an initial decline m ozone levels due to reaction number 8 . This 
decline is an indication that an initial ozone concentration of 0.050 ppm is 
incompatible with initial N0 X in the 0.020 - 0.050 ppm range and with 75 
percent of this N0 X being NO. Calculations with initial NO 2 taken to be 80 
percent of the initial N0 X showed no initial decline for any simulation. 

The results of calculations with this ratio are presented m Appendices A 
and B, with and without mixed layer rise. 

The lowest NMHC case, figure 2a, indicates that 0.100 ppmC does not 
represent sufficient NMHC to discriminate between various N0 X levels; in 
fact, an unexpected inverse trend between O 3 and initial N0 X is displayed. 

No doubt, part of the explanation of the inverse trend lies in the 
establishment of a consistent O 3 /NO/NO 2 pattern, that is, a mixture that 
would not undergo a rapid change. The 0.200 ppmC NMHC case, figure 2b, 
displays a direct relationship between O 3 and N0 X after 2 hours of 
simulation, except for the highest N0 X levels. Figure 2 (c and d) , 0.300 
and 0.400 ppmC of NMHC, displays a direct and increasingly severe 
relationship between O 3 and initial N0 X . With NMHC greater than 0.200 
ppmC, N0 X m the range 0.020 - 0.050 ppm produces the expectation of an 
afternoon oxidant episode in the context of the model. 

Figure 3 (a, b, c, d) is a cross-plotting of what was presented m 
figure 2. The different plots show the change in ozone concentration with 
time and varying initial NMHC concentration, with a different initial N0 X 
concentration for each plot. Except for the lowest N0 X concentration, 
figure 3a, and highest NMHC concentration (and then only after about 500 
minutes) for this figure, there is a direct relationship between O 3 and 
initial NMHC concentration. An especially severe relationship exists at 
the highest N0 X level, figure 3d ( 0.050 ppm N0 X ) . 

While an inverse relationship can exist between O 3 production and 
initial NMHC and N0 X concentration, the predominant relationship for the 
concentration range of the simulation (table I) is a direct one. One 
could also inquire if O 3 production is related to the NMHC/N0 x ratio and 
could observe intable- 1 three cases in which this ratio is 10.0. The 
corresponding curves, however, in figures 2 or 3 show no particular 
significance for the ratio of NMHC to N0 X within the range of concentrations 
used in the simulation. 


Traveling Air Parcel Model 


Using an intermediate case as a base, sources were added during the 
first 2 hours of simulation as a representation of an air parcel passing 
over an urban area. The intermediate case consisted of the following ini- 
tial conditions: NG X 0.030 ppm, NMHC 0.200 ppmC, and O 3 0.050 ppm. These 

concentrations are probably high compared to normal background levels. 
Figure 4 displays the variation of ozone concentration with time for a 
NMHC loading of 108.8 yg m - ^ hr - "* - * and with the N0^ loading varying from 


* 

1 . 0 ppm 


( 1 . 2187 x 1 0 1 *) 


(mole w eight ) ^pg 


(T,K0 
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5.UU to 21.8 yg nT^ hr - "* - . Figure 5 displays the variation of ozone 
concentration with time for a NOx loading of 5.44 yg m - ^ hr - "* - and with the 
NMHC loading varying from 54.4 to 217.6 yg m - 3 hr - -*-. A NMHC loading at the 
54.4 yg m“3 hr - -*- level is consistent with the yearly regional loading of 
100,000 tons. The N0 X loading was varied over a range of values to produce 
a change in ozone production. Any increase in N0 X or NMHC over the initial 
range of concentrations produces an increase m photochemical oxidant over 
the previously calculated case with the exception of the highest NMHC 
loadings of figure 5. 


CONCLUDING REMARKS 

Computer predictions of ozone concentration for 10-hour photochemical 
simulations, beginning at 0800 l.d.t., have been calculated as a function of 
initial N0 X and nonmethane hydrocarbon (NMHC) input concentrations. A 
contemporary chemical mechanism (including diurnally varied rate 
coefficients for six photochemical decompositions) and procedure (which 
specifies, in part, an assumed constitution of N0 X and NMHC) were used in 
these photochemical box-model calculations. The range of precursor 
species concentrations used m the simulations (N0 X varying from 0.020 to 
0.050 ppm and NMHC varying from 0.100 to 0.400 ppmC) is generally thought 
to span that characteristic of the Southeastern Virginia region on a summer 
morning before an afternoon photochemical oxidant episode. With the 
exception of the lowest values for the NMHC, afternoon photochemical 
oxidant concentrations greater than 0.010 ppm are predicted by the model. 

All photochemical box-model calculations predict a rising ozone level from 
morning to afternoon with the concentration varying between 0.040 and 0.200 
ppm, which is consistent with observations m the Southeastern Virginia 
region. It should be noted, however, that the present calculations present 
a worst case (highest ozone) situation since the contemporary procedure 
assumes a lower dilution rate than is probable. That the chamber produced 
chemical scheme yields predictions consistent with observations for such a 
crude model lends encouragement to the development of models with more 
sophisticated meteorology and to the obtaining of more exact species 
concentration data. 

Since the initial input to the model, which begins m the morning with 
rural air, is necessarily aged air with pollutant precursors generated the 
previous day, it is thought that the initial N0 X should be characterized by 
a higher proportion of NO 2 than that used m the simulation (NO 2 • NO = 1:3). 
The computer simulation is consistent with this idea as is shown by the 
initial decline in ozone for all cases. However, peak ozone is not affected 
significantly by the initial proportions of N0 X . As a corollary thought, 
since n-butane is less reactive with oxidants m air than propylene, the 
proportion of propylene to n-butane, representing the initial input of NMHC 
to the model, should be lowered such that there is a greater proportion of 
n-butane in the representation of NMHC. This would result m lowered ozone 
production since the revised NMHC mixture would be less reactive. 

Simulations with the Lagrangian (traveling air parcel) model yield 
relatively higher ozone levels than those with photochemical box model 
since they are, m effect, a coupling of an intermediate box-model 
calculation with 2-hour source additions of N0 X and NMHC. The mid-range of 
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the assumed initial precursor concentrations has such a potential for oxidant 
production that any additional contribution of either N0 X or NMHC would lead 
to significant increases in ozone. 

The initial N0 X and NMHC concentrations were represented in the 
calculations by certain ratios of NO2 to NO and propylene to n-butane, 
respectively, as were the source additions (though not necessarily by the 
same ratios). It is felt that the use of these assumed ratios can represent 
a significant source of error, with the use justified only by its allowing 
calculation with insufficient data; for example, hydrocarbon data are 
available only as methane and total hydrocarbon concentrations, not 
concentrations of individual hydrocarbons. In view of the development of 
detailed chemical mechanisms for the hydrocarbons, it is believed, along 
with others (see, for example, A. C. Lloyd, NASA RP 1022, pp. 278-279, 

1978) (ref. 13), that the requirement exists for much more detailed air 
quality data measurements. It is also clearly evident that other critical 
inputs are accurate solar insolation, emissions inventory, and meteorological 
dat a . 
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APPENDIX A 


The changes in ozone concentration with time for various initial N0 X 
concentrations are shown in this section. These calculations were made 
assuming a N 02 /N 0 X ratio of 0.80. This ratio of NO 2 to N0 X did not 
produce the initial decline of ozone during the first few minutes of the 
simulation. The implication is that the original NO 2 to N0 X ratio used 
was not consistent with the initial assumption of 0.050 ppm O 3 . 
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Figure Al. - Change of 0^ concentration with time for various initial 
NO concentrations. 
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Figure A3. - Change of 0 3 concentration with time for various initial 
NO concentrations. 
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APPENDIX B 


The changes in ozone concentration with time for various initial N0 X 
concentrations are shown in this section. These calculations were made 
assuming a N02/N0 X ratio of 0.80; however, dilution from mixed layer 
height rise was not included. With no dilution, the levels of ozone 
concentration are higher than cases where the mixing layer height was 
allowed to change. 


15 



Ozone, ppm 



Figure Bl. - Change of 0^ concentration with time for various initial 
NO x concentrations. 


16 



Ozone, ppm 



Time, Local Daylight Time 


Figure B2. - Change of 0^ concentration with time for various initial 
NO w concentrations. 
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Figure B3. - Change of O^ concentration with time for various initial 
NO concentrations. 

X 


18 




Ozone, ppm 


0.25 


0.20H 


0.15b- 


0.10 b- 


Q 05 


INITIAL NMHC 0.400 ppmC 
NO z /NO x =0.80 
NO DILUTION 


| t , | iii '•>• m m m m m.M 



INITIAL NO. 

o 0. 020 ppm 
□ 0. 030 ppm 
O 0. 040 ppm 
A 0. 050 ppm 


i i i i i i i i i i i 


i I i i i I i i 


i i I i i i I i i i 


0800 


1000 


1200 1400 

Time, Local Daylight Time 


1600 


1800 


Figure B4. - Change of 0^ concentration with time for various initial 

NO concentrations, 
x 

19 


H 




REFERENCES 


1. EPA-Released State Air Quality Data Show Areas Not Meeting Federal 

Standards. Journal of the Air Pollution Control Association, vol. 28, 
no. U, Apr. 1978, pp. 378-380. 

2. DeMarrais, Gerard A.: The Ozone Problem m the Norfolk, Virginia Area. 

EPA-600A-78-006, Jan. 1978. 

3. Wagner, H. S., Gregory, G. L. , and Buglia, J. J. : The Southeastern 

Virginia Urban Plume, A Test Site for Remote Sensors. Seventy-first 
Annual Meeting of the Air Pollution Control Association, Houston, TX, 
June 25-29, 1978. 

k. Gear, C. William: Numerical Initial Value Problems in Ordinary 

Differential Equations. Prentice-Hall, Inc., cl971. 

5. Hecht, Thomas A.; Seinfeld, John H. , and Dodge, Marcia C.. Further 

Development of Generalized Kinetic Mechanism for Photochemical Smog. 
Environmental Science and Technology, vol. 8, no. 4, Apr. 197*+, 

pp. 327-339- 

6. Dodge, M. C. and Hecht, T. A. Rate Constant Measurements Needed to 

Improve a General Kinetic Mechanism for Photochemical Smog. 
International Journal of Chemical Kinetics, Symposium No. 1, 1975. 

7. Schere, K. L. and Demerjian, K. L. A Photochemical Box Model for Urban 

Air Quality Simulation. Proceedings of the Fourth Joint Conference on 
Sensing of Environmental Pollutants, New Orleans, LA, Nov. 6-11, 1977, 
pp. 1+27-H33. 

8. Meyers, E. L. ,Jr. ; Summerhays, J.E. ; and Freas, W. P. • Uses, Limitations 

and Technical Basis of Procedures for Quantifying Relationships 
Between Photochemical Oxidants and Precursors. U.S. EPA Publication 
No. EPA-U 50 / 2-77- 021a, Nov. 1977- 

9. Schere, Kenneth, L. and Demerjian, Kenneth L. Calculation of Selected 

Photolytic Rate Constants over a Diurnal Range. EPA-600 /*i-77-015 , 

Mar. 1977- 

10. Demerjian, K. L. . Photochemical Air Quality Simulation Modeling: 

Current Status and Future Prospects. International Conference on 
Photochemical Oxidant Pollution and Its Control Proceedings, Raleigh, 
NC, vol. I, Sept. 12-17, 1976. 

11. Graedel, T. E. , Farrow, L. A., and Weber, T. A. Kinetic Studies of 

the Photochemistry of the Urban Troposphere. Atmospheric Environment, 
vol. 10, 1976, pp. 1095-1116. 

12. Dodge, Marcia C. Effect of Selected Parameters on Predictions of a 

Photochemical Model. EPA-600/3-77-0*+8 , June 1977. 

13. Lloyd, Alan C. The Role of Computer Modeling of Photochemical Smog in 

Defining Existing Measurement Needs. NASA RP-1022, 1978, pp. 278-279- 


20 






Ozone, ppm 


INITIAL NMHC 0.100 ppmC 


i i i I i i i I i i i 1 i i i 


INITIAL NO 
) 

o 0.020 ppm 
□ 0.030 ppm 
O 0.040 ppm 
A 0.050 ppm 

i i I i i i I i 


Time, Local Daylight Time 

Figure 2a.- Change of O^ concentration with time for various initial 
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Figure 2c.- Change of 0^ concentration with time for various initial 
NO concentrations. 
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Figure 3a.- Change of 0^ concentration with time for various initial 
NMHC concentrations. 
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Figure 3b.- Change of 0^ concentration with time for various initial 
NMHC concentrations. 
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Figure 3c.- Change of 0^ concentration with time for various initial 
NMHC concentrations. 
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